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I.  INTRODUCTION 


A.  PROBLEM  DESCRIPTION 

It  is  common  in  the  real  world  to  use  a  motor  to  drive  a 
mechanical  load.  So,  waen  the  structure  is  subjected  to 
periodic  forcing,  it  may  well  exhibit  a  resonance  or  a 
number  of  resonances  (harmonics)  wnich  can  lead  the  system 
to  instability.  Two  typical  examples  of  the  above  problem 
are  the  head  servo  in  a  disk  memory  and  the  electric  typing 
machine  used  with  word  processors. 

In  order  to  be  familiar  with  tie  problem  and  see  what 
causes  our  system  to  instability  we  have  to  study  the  time 
and  frequency  response  of  a  system  which  becomes  unstable 
due  to  mechanical  resonances. 

In  the  following  paragraph  we  consider  the  problem  in 
the  frequency  domain  where  it  is  easier  to  define  and  under¬ 
stand. 

B.  FREQUENCY  DOMAIN 

Considering  a  typical  control  system  we  know  that  in 
order  for  this  system  to  be  stable  we  need  a  positive  phase 
margin.  The  positive  phase  margin  happens  if  the  gain  cross¬ 
over  of  the  system  is  at  a  lower  frequency  than  the  phase 


crossover  or.  the  3ode  diagram  is  shown  in  Figure  1.1.  So  one 
resulting  phase  margin  vector  ('y)  is  positive. 

If  now  we  attempt  to  drive  one  or  more  mechanical  loads 
with  our  control  system  than  it  is  possible  that  the  system 
may  become  unstable  due  to  the  mechanical  loads.  The  reason 
for  this  instability  is  tie  rssonaaoe  peak  or  peaks  created 
by  the  mechanical  loads.  In  general  we  can  summarise  tha 
problem  by  the  following  two  cases. 

(1)  At  the  resonances  tha  peak  of  the  gain  curve  does  not 
rise  above  the  zero  d3  axis. 

(2)  at  the  resonanses  tha  pear  of  the  gain  risas  above  tha 
zero  dB  axis 

Analyzing  each  case  separately,  wa  can  easily  understand 
that  for  the  first  case  chare  is  no  stability  problem  since 
tha  resonances  do  not  excaad  the  zero  dB  axis,  which  means 
that  the  phase  margin  we  had  remains  the  same  and  if  our 
system  was  stable  it  remains  stable. 

Considering  now  the  second  casa  wa  see  that  if  one  of 
tha  resonance  peaks  exceeds  the  zero  dB  axis  it  will  create 
a  new  gain  crossover  which  can  make  the  system  unstable  if 
tha  resulting  phase  crossover  is  at  a  lower  frequency  than 
tha  gain  crossover  i.e.  tne  resulting  phase  margin  will  be 
negative. 


1. 1 


Possible  Phase  M 


Assuming  that  the  reader  now  is  familiar  with  the 
mechanical  resonances  problem  we  cm  proceed  ar.d  try  to 
solve  the  problem  i.e.  to  compensate  a  system  which  has  bean 
unstable  due  to  these  resonanses,  using  two  different 
met  hods. 

C.  THE  HOD EL 

The  model  which  we  will  use  in  the  major  part  of  the 
srudy  has  to  be  initially  stable.  Sa  we  choose  conveniently 
a  third  degree  system  with  one  zero  and  three  poles  with 
locations  that  night  appear  in  a  physical  system  whose 
Laplace  transform  has  the  following  form: 

Qrsl_2500(,Q.ISH')  (i3=  1.1, 

S(0.2S+ 0(00033  5+ 1) 

As  we  see  in  Figures  1.2  ,  1.3  the  frequency  response  of 
the  system  proves  that  the  system  is  stable  since  we  have  a 
positive  phase  margin  of  '42  degrees.  In  Figure  1.4  we  can 
sea  the  time  response  of  the  system  with  a  small  oscillation 
for  0.04  seconds,  which  also  proves  the  stability  of  the 
system. 
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Figure  1.2  Gain  Oarve  of  tha  Stable  Systee 


Figure  1.3  Phase  Curve  of  tie  stable  Model 


OUTPUT 


Now  we  have  to  add  two  typical  resonances  in.  or  far  to 
sea  how  the  system  will  become  unstable. 

The  two  selected  mechanical  resonances  are  as  follows: 
FIEST  HAfiHONIC  at  frequency  f1=2000  HZ  with  lamping 
ratio  £*=0.005 


G  (s) 


_ l 

0.25 105'+ 0.510% +1 


(Egn  1.2) 


SECOND  HAfiHONIC  at  frequency  £2=5000  HZ  wirr.  camping 
ratio  £”=0.0000  158 


G($)= - 5 - - - i -  (Sgn  1.3) 

0.A  I0' 5*4  0.^10  5  +  1 

After  the  introduction  of  the  two  mechanical  resonanses 
the  transfer  function  becomes: 


G(s)= 


2500 


Co.lS+f) 


sfo.2S4  (Ko.035StiXo.25-  I0'V+  0.5-  lO'%+^Xo.4  IOS+a63 10%+i) 


(Eqn  1.4) 


The  following  pages  present  frequency  responses  on 
Figures  1.5,  1.6  and  time  response  on  Figure  1.7  of  this 

system,  including  the  two  resonances.  . 

From  this  point  in  this  Thesis  we  will  call  Gu  (s) 

•  Model' wherever  it  is  needed.  The  '  u»  means  uncompensated. 

In  Appendix  the  programs  mostly  used  in  this  Thesis  ace 
given.  Program  1  is  suitable  for  culcalating  the  frequency 
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response  of  the  system  with  suitable  modif ication  for  each 
case.  Frogram  2  is  suitable  for  calculating  the  time 
response  of  the  system  with  suitable  modification  for  each 
case  also.  Program  3  is  suitable  for  calculating  the  time 
response  of  the  system,  whan  wa  usa  is  a  compensation  method 
imaginary  zeros  and  complex  poles.  The  graphs  have  been 
obtained  using  the  Vacsatec  plotter  of  the  Naval 
Postgraduate  School  computer  IBM-3333  using  the  Digital 
Simulation  Package  (DSL)  and/or  S3 3?  III.  In  the  phase 
curves  when  the  phase  exceeds  the  —353  degress  the  plotter 
jumps  tc  its  next  value  aid  the  curve  looks  discontinued  but 
in  reality  it  is  not. 
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Pigure  1.6  Phase  Curve  of  the  Systen  with  the  Resonances 


OUTPUT 
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Figure  1.7  Tiae  Response  of  the  systea  with  the  Resonances 


D.  DEFINITIONS 


In  order  to  follow  the  thesis  easily,  it  is  a  good  idea 
to  refresh  our  memory  with  a  few  definitions. 

1.  GAIN  CROSSOVER  :  This  is  the  point  on  the  plot  of  the 
transfer  function  at  which  the  magnitude  is  unity  i.e.  the 
point  at  which  the  gain  curve  crosses  the  zero  iB  axis. 

2.  PHASE  CROSSOVER  :  This  is  the  point  on  the  plot  of 
the  transfer  function  at  which  the  phase  is  -180  decrees 
i.e.  the  point  at  which  the  phase  curve  crosses  the  -180 
degrees  axis. 

3.  PHASE  MARGIN  :  This  is  defined  by  the  equation: 

~Y  =180+|G(s)|  and  detscnines  the  amount  of  phase  shift 
required  to  place  the  system  at  the  stabiity  limit. 

4.  ROOT  LOCOS  :  Is  a  plot  of  the  roots  of  the  character¬ 
istic  equation  of  the  closed  loop  system  as  a  function  of 
the  gain  or  some  other  variable  paraneter. 

5.  RESONANCE  PEAK  :  Is  the  maximum  amplitude  of  the  gain 
curve.  Of  the  frequency  response.  It  occurs  at  the  resonant 
frequency,  Wr,  and  is  due  to  complex  poles  (roots). 

6.  RESONANT  FREQUENT?  (Nr):  The  frequency  at  which 
|G(jw)|  has  its  peak  value 

7.  DAMPING  RATIO  :  Is  the  factor  £“=F/Fc  where  F  is  the 
actual  damping  and  Fc  =  2 rVj xtc  the  critical  factor  of  a 
quadratic  having  form:  G(s)«  K/(Jsa+FS+lO 


II.  COMPENS&II CH  USING  IS  AGINARY  ZEROS  AND  COMPLEX  POLES 

A.  HETHOD  DESCRIPTION 

In  this  method  we  will  use  a  compensator  which  includes 
imaginary  zeros  and  complex  poles.  This  compensator  has  the 
following  general  form: 

( J  (5)  = - — -  (2gn  2.1) 

5a+ 

Where  a  is  the  frequency  of  the  puce  zero  on  the  imaginary 
axis  and  Wn  is  the  frequency  of  the  complex  poles. 

The  philosophy  of  cats  compensator  is  that  we  cry  to 
raise  the  phase  curve  at  a  frequency  slightly  lower  than  the 
resonant  frequency  in  order  to  avoid  the  negative  phase 
margin.  Then  we  lower  the  phase  curve  after  the  resonance 
has  occuraa  so  that  the  gain  curve  lies  under  this  phase 
peak  and  the  resulting  new  phase  margins  are  positive  at 
both  gain  crossovers,  in  other  words  we  try  to  reshape  the 
phase  curve.  For  this  reason  we  put  the  zero  in  the  vicinity 
of  the  resonant  frequency  (1.030  Hz)  and  specifically  at  a 
slightly  lower  frequency  (1  950  Hz).  In  this  way  we  try  to 
raise  the  phase  curve  before  the  resonance  happens.  Then  in 
order  to  return  the  system  to  its  original  condition  we 
introduce  the  complex  poles  at  a  frequency  of  2350  Hz. 
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In  the  following  pages  we  see  the  frequency  response  and 
the  time  response  of  the  'Sods!*  when  we  apply  this  method. 

B.  METHOD  APPLICAI  ION 

The  first  data  we  use  for  our  compensator  is  a  frequency 
1950  Hz  for  the  imaginary  zero  and  a  slightly  greater 
frequency  for  the  complex  poles,  2050  Hz.  The  resonant 
frequency  is  2000  Hz.  Keeping  the  above  data  constant  we 
vary  the  damping  ratio  of  the  compensator  in  order  to  see 
how  this  variation  affects  the  filter  . 

In  Figures  2.1,  2.2,  2.3  we  lavs  the  gain  and  phase 

curves  and  the  Time  response  of  the  system  wirh  damping 
ratio  £*=0.05.  In  Figures  2.4,  2.5,  2.6  we  have  the  gain, 

phase  and  time  response  of  the  system  with  damping  ratio 

1=0.  1. 
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Figure  2.2  Phase  :urv=  with  ^=0.05. 
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Pigure  2.5  Phase  Curve  with  £"*0.1 


Observing  Figures  2.1  through  2.5 


:  v  i  o  a  s 


method  has  not  worked  since  the  system  remains  unstable.  The 
two  different  damping  ratios  used  i.e.  0.05  and  0.1  are  the 
two  representatives  for  this  point  since  we  do  not  want  to 
increase  the  size  of  the  thesis  with  unimportant  graphs,  but 
in  reality  we  tried  many  more,  concluding  that  the  damping 
ratio  does  not  have  much  effect  on  the  system  at  this  point 
of  the  study. 

If  we  observe  the  graphs  closer  or  we  overlap  the  gain 
curves  over  the  phase  curves  we  will  see  that  we  have  three 
gain  crossovers  and  consequently  three  phase  margins.  The 
first  is  the  phase  margin  of  the  original  model  before  the 
resonances  and  it  is  stable.  The  other  two  are  created  from 
the  gain  crossover  of  the  first  peak  of  the  •  Siodel*  which 
crosses  the  zero  d3  axis  at  two  points.  The  first  of  these 
two  phase  margins  is  positive  and  the  other  is  negative. 
That  means  something  happened  because  before  we  used  the 
filter  both  phase  margins  were  negative,  so  the  filter  has 
worked  but  partialy  and  specifically  the  frequency  of  the 
zeros  was  correct  since  the  phase  curve  has  been  raised  and 
the  second  of  the  phase  margins  became  positive.  What  we 
have  to  do  now  is  to  adjust  the  frequency  of  the  poles  in 
order  to  reshape  the  high  frequeny  part  of  the  phase  curve 
and  have  the  third  phase  margin  positive  too,  which  means 
that  the  system  will  be  stable.  In  Figures  2.7  through  2.15 
we  can  see  the  gain,  phase  and  time  response  of  the  system 


at  pole  frequencies  2400, 2600, 2800  Hz  respectively.  t^3 
frequency  of  the  zeros  is  also  reduce!  at  1900  Hz  i.e.  50  Hz 
less  in  order  to  see  if  we  can  have  a  small  variation  of  the 
zeros  for  application  purposes. 
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Pigure  2.8  Phase  Curwa  with  Poia  at  Frequency  2400  Hz. 
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Figure  2.14  phase  "arve  with  poles  at  2890  Hz  . 
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After  examination  of  Figures  2.7  through  2.15  we  observe 
that  the  system  becomes  stable  for  pile  frequencies  2500  and 
2830  Hz  respectively  with  the  only  difference  being  that  the 
system  has  less  transient  oscillation  when  we  put  the  poles 
at  frequency  2800  Hz  (system  oscillate  oscillates  for  0.25 
seconds)  than  at  frequency  2500  Hz  (system  oscillates  for 
0.5  seconds) 

In  order  to  understand  more  completely  the  behavior  of 
the  compensated  system  for  various  gains  it  is  wrse  to 
obtain  the  root  locus  of  tae  system.  3o  we  have  to  calculate 
the  characteristic  equation  of  the  system.  The  character¬ 
istic  equation  comes  from  the  transfer  function  of  the 
system  with  the  following  formula  3(s)*1=0  but 
G  (s )  =N  (s) /D  (s)  where  N(s)  is  the  numerator  and  D(s)  is  the 
denominator  of  the  transfer  function  and  finally  the 
CHARACTERISTIC  EQUATION  =N(s)+D(s). 

After  the  calculations  we  obtain  the  following  equation 
of  ninth  degree: 

C.E  =  S+885.I58S8+ 37.02  loVf  2.8l3-ldV+  3.32-10 S5+ 

+  I.59-I0V+  S3(  8.028' 10*°+  5.257.  fOls«)+  £(2.39  ■  10^4 

4  3.257. 10,7k:)+ S(U7G- left.  1.176- lo”k)  4  (Eqn  2*2) 

+  I.  I7G  •  I07*k. 

In  Figures  2.16  we  have  plotted  the  roots  of  the  system 
varying  the  gain  (K)  from  3.01  to  333444  . 
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From  the  root  locus  of  the  system  we  observe  *  a-  t.-.e 
system  can  have  roots  in  the  right  naif  plane,  meaning  that 
for  some  specific  values  of  gain  ( K)  the  system  becomes 
unstable. 

Summarizing  our  results  up  to  this  point  we  can  say  that 
this  method  will  work  succesfully  and  generally  we  can 
stabilize  a  system  including  mechanical  resonances  with  the 
proposed  filter  putting  imaginary  zeros  in  the  vicinity  of 
the  resonant  frequency  ana  specifically  100-150  Hz  below  tne 
frequency  of  the  system  complex  poles.  The  complex  poles  of 
the  filter  are  placed  a  few  hundred  Hertz  above  the  resonant 
frequency.  Choice  of  the  filter  poles  depends  also  on  the 
frequencies  of  any  additional  resonances  in  the  system,  and 
trial  and  error  methods  are  needed  to  find  the  best  pole 
location.  It  is  wise  aiso  to  focus  our  attention  or.  the 
system  gain  (KJ  when  we  apply  this  method  in  order  to  avoid 
transition  of  other  roots  into  the  right  half  plane. 
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III.  COMPENSATION  CJSIN3  I3MPLEX  POLES 


A.  METHOD  DESCRIPTION 

When  a  system  is  unstable  because  of  a  resonance,  the 
bcie  diagram  shows  that  the  peaking  of  the  magnitude  curve 
generates  two  additional  gain  crossover  points,  while  the 
rapid  decrease  in  phase  due  to  the  complex  poles  causes  The 
phase  curve  to  pass  through  an  odd  multiple  of  7T  (-180  or 
-540  degrees)  at  a  frequency  between  these  gain  crossover 
points.  Thus,  at  the  higher  frequency  gain  crossover  there 
is  a  negative  phase  margin  indicating  that  the  system  is 
unstable.  In  using  an  additional  pair  of  complex  poles  to 
stabilize  the  system,  we  reshape  the  pnase  curve  so  that  the 
phase  crossover  occurs  at  a  frequency  lower  than  that  of 
either  gain  crossover  due  to  the  resonance.  This  eliminates 
the  possibility  of  the  negative  phase  margin  (i.  e.  elimi¬ 
nates  the  encirclement  of  the  Nyquist  point)  .  Note  that  in 
designing  the  compensator  the  reshaping  of  the  gain  curve 
(due  tc  the  compensator  complex  poles)  must  not  generate  any 
additional  gain  crossover  points.  This  means  that  the 
damping  ratio,  ,  of  the  compensator  poles  must  not  be  too 
small.  Some  compromise  may  therefor  be  expected,  since  we 


would  prefer  a  very  sharp  decrease  in  phase  which  in  turn 
requires  small  £*• 


B.  METHOD  APPLICATION 


It  is  desirable  to  live  c  he  plots  for  the  system  in 
order  to  realise  how  the  filter  works.  In  the  following 
Figures  we  have  the  Nyguist  plots  of  the  stable  system 
without  the  resonances  and  with  the  resonances.  Figure  3.1 
shows  the  Nyquist  plot  of  the  system  before  the  introduction 
of  the  resonances.  This  plot  is  in  large  scale  in  order  to 
show  ail  the  values  of  tie  imaginary  versus  real  coordi¬ 
nates,  but  it  is  not  clear  whether  tie  -1  point  is  inside  or 
outside  the  curve.  In  Figure  3.2  we  can  see  than  the  -1 
point  is  not  included  in  the  curve.  In  Figure  3.3  we  can  see 
the  system  after  the  introduction  of  resonances.  If  is  easy 
to  observe  that  the  -1  point  is  inside  the  curve  which  makes 
the  system  unstable. 


Figure  3.3  The  Systao  with  the  Resonance  Peaks. 


Now  we  will  try 


*he  system 


to  s*aaiiise 

ccaplex  pole  filter.  The  filter  will  have  tie  following 
transfer  function: 

- - -  (Eqn  3.1) 

5’+  2fcj„s  +  GJ,a 

Where  we  have  to  specify  the  pole  location  ana  the  damping 
ratio.  From  the  experience  of  the  previous  chapter  we  .know 
that  =0.1  as  a  reasonaoie  value  waini  does  not  create  addi¬ 
tional  peaks.  So  we  will  use  this  vaLae  and  we  will  vary  the 
pole  location  by  trial  and  error. 

Since  we  want  the  phase  lag  to 'be  increased  to  more  than 
-360  degrees  what  we  think  first  is  to  locate  the  poles  in 
the  vicinity  of  the  first  resonance  and  see  if  this  choice 
stabilizes  the  system.  If  that  happens  we  will  change  the 
pole  location  to  lower  and  higher  frequencies  from  the  reso¬ 
nance  peak  in  order  to  see  how  flexiole  the  method  is,  since 
in  the  next  chapter  we  are  planning  to  compere  the  two 
met  hods. 

As  a  first  approach  we  will  use  pole  frequency  1950  Hz 
i.a.  50  Hz  lower  than  the  resonance  peak  with  £"=0.1  .  In 


the  following  Figures  we  can  see  the  gain,  the  phase  and  th° 
time  response  of  the  system  using  this  compensator. 


XSCRLE=  0.20  UNITS/INCH  FREQ  (H7>)  BUN  NO.  1 

TSCflLE*  140.00  UNITS/INCH  '  PLOT  NO.  2 


Figure  3-5  Phase  Curve  oE  the  System  with  Poles  at  1950  Hz 


XSCfllE*  O.0'4  UNITS/INCH  Tmr  ferr>  RUN  NO.  1 

YSCfflE=  0.40  UNITS/INCH  '  PLOT  NO.  1 


Figure  3.6  Ti*e  aesponce  of  the  System  with  Poles  at  1950  Sr 


Observing  Figures  3.4  through  3.3  we  see  that  the  filter 
has  worked  as  expected.  The  phase  lag  has  exceeded  -360 
degrees.  From  the  time  response  curve  of  Figure  3.6  we  can 
see  a  small  oscillation  for  about  0.  OS  seconds  and  then  the 
system  becomes  stable.  The  Nyquist  plots  of  Figures  3.7  and 
3.3  verify  the  stability  of  the  system  too.  In  Figure  3.3  we 
can  see  the  Uyquist  plot  magnified  specifically  at  the 
vicinity  of  the  -1  point  and  it  is  easy  to  observe  that  the 
point  is  lied  outside  tae  curve. 

After  we  have  proved  that  tha  method  has  worked  sucess- 
fuly  it  is  necessary  to  find  the  limits  at  the  frequency 
axis  we  can  locate  the  complex  poles  in  order  for  the  system 
to  become  stable.  First  we  will  remove  the  poles  to  higher 
frequencies  and  then  to  lower  frequencies. 

Figures  appear  the  results  of  that  research. 


In  the  following 


cea) 


XSCflLE=  0.20  UNITS/INCH  FREQ  (HZ)  RUN  NO.  1 

YSCALE=  8.00  UNITS/INCH  '  '  PLOT  NO.  1 


Figure  3.9  Gain  Curve  ot  the  Systaa  with  Poles  at  1900  Hz 
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TRANSIENT  RESPONSE  OF  THE  COMPENSATED  SYSTEM 
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Figure  3.11  Tiae  Eesponss  of  the  Srstea  with  Poles  at  1900  Hz 


xscale*  0.08  uNiTs/rNCH  FREQ  (Hz)  RUN  N0-  1 

YSCflLE*  8.00  UNITS/INCH  PLOT  NO.  1 


Figure  3.12  Gain  Curve  of  the  Systai  with  Poles  at  2300  Hz 
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14  Tiae  Besponsa  of  the  System  with  Poles  at  2300  Hz. 
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Pigure  3.18  Gain  Cuevas  for  Paribas  Pole  Freguencies. 
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In  Figures  3.3  through  3.11  we  ran  see  what  barrens  if 
we  reduce  the  pole  location  53  Hz.  3o  we  put  the  poles  at 
1900  Hz  and  what  we  see  is  that  the  system  becomes  unstable 
as  shown  in  Figure  3.11#  the  time  response  of  the  system. 
That  indicates  that  the  lower  limit  of  the  pole  location  is 
at  a  frequency  about  100  Hz  lower  than  the  resonant 
frequency.  Now  we  try  to  find  the  upper  frequency  limit,  in 
Figures  3.12  through  3.14  we  can  see  that  the  system  is 
stable  as  we  move  the  poles  to  higher  frequencies. 
Specifically  if  we  observe  Figure  3.13  we  see  that  there  is 
a  longer  oscillation  period  (0.24  seconds)  compared  with  the 
oscillation  period  of  the  stabilised  system  with  poles 
located  at  1950  Hz  which  was  0.09  seconds.  This  is  one 
indication  that  when  we  increase  the  frequency  of  the  poles 
we  have  a  larger  oscillation  time.  In  Figures  3.15  thzougn 
3.17  we  place  the  poles  at  2730  Hz  and  the  oscillation  time 
has  further  increased  to  3.53  seconds  which  certifies  again 
what  we  claimed  above. 

In  the  last  two  Figures  we  can  see  the  change  in  shape 
of  the  gain  and  phase  graphs  as  we  increase  the  frequency  of 
the  pole  location  to  2200,  2300,  2303  Hz,  respectively. 

Closinq  this  chapter  we  conclude  that  this  method  worked 
for  our  'Model*.  Specifically  we  can  say  that  when  we  locate 
the  poles  of  the  filter  we  must  take  into  acouat  the  oscil¬ 
lation  period  which  is  /ery  important  for  a  system.  The 
range  of  frequency  with  which  we  cm  locate  the  poles  is 


3 


IV.  C0HCLPSI3NS  AND  8 SC D S EN D A TIONS 


A.  CONCLUSIONS 

Two  methods  of  compensation,  for  instability  due  to  reso¬ 
nances  have  been  studied.  Me  can  conclude  that  both  were 
sucessful  and  both  can  be  used  to  compensate  systems 
including  mechanical  resonances.  In  the  first  netnod,  using 
pure  imaginary  zeros  and  complex  poles  we  need  to  locate  the 
zeros  at  a  frequency  100-  150  Hz  lower  than  the  resonnce 
frequency  of  the  system,  and  the  complex  poles  a  few  hundred 
Hz  higher.  This  method  has  the  disadvantage  that  it  gives  an 
oscillation  period  that  is  undesirably  long  and  additionally 
there  are  some  gain  values  that  which  make  the  system 
unstable.  The  second  method  using  only  complex  poles  appears 
more  succesful  because  we  have  a  larger  range  of  permissible 
pole  location.  The  oscillation  period  can  be  signficantly 
reduced  if  we  build  the  filter  poles  with  frequency  near  the 
resonant  frequency. 

B.  8  ECOHHENDATIOHS 

In  this  study  we  used  only  a  ' Sodel*  built  for  this 
problem.  Before  we  conclude  that  tie  methods  are  useful  in 
applications  it  is  wise  to  try  to  apply  both  methods  to  more 
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APPENDIX  k 
COM  PCI  I  2R  PROGRAMS 

This  appendix  is  composed  of  listings  of  the  various 
DSL/360  and  CSMP  III  programs  used  in  •'■he  computer  simula¬ 
tions  throughout  this  study. 

Specifically  on  page  7 3  is  the  program  used  to  obtain 
the  frequency  response  of  tha  'Model'  without  the  reso¬ 
nances.  On  page  74  is  the  program  used  to  obtain  the 
frequency  responce  of  the  'Model'  with  tha  resonances.  On 
page  75  is  the  program  used  to  obtain  the  frequency  response 
of  the  compensated  system  using  complex  poles  and  imaginary 
zecos.  On  page  76 -is  tha  program  used  to  obtain  the  time 
responseof  the  compensated  system  using  complex  poles  and 
imaginary  zeros.  On  page  77  is  the  program  used  to  obtain 
tha  frequency  response  of  the  compensated  system  using  only 
complex  poles.  On  page  78  is  the  program  used  to  obtain  the 
time  response  of  the  compensated  system  using  only  complex 
poles.  To  obtain  the  root  locus  of  page  43  we  used  the 
subroutine  HOOTLO  on  the  CBM  3  033  of  the  Naval  Postgraduate 
School.  On  page  79  is  tha  program  used  to  obtain  Figures 
3.18  and  3.19. 
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T  ITl  5  FREQ'JEN  C  V  RESPONSE 
TITLE  MODEL  WITHGUT  RESONANCES 
D  COMPLEX  S,N,D»G 

D  COMMON/ CARE A/ S  »N, C »G 

INT3ER  N  PLOT 
CONST  NPLQT=  1 

CONTRL  FINTIM=4.0t  CELT=C. C0C6»DELS=0. C025 
PRINT  4.0B-3 ,W,PHI »RE, IM,MD8,MAG 
DYNAMIC 

RW=RAMP ( 0.0 ) 

LCGW=RW 
W=10.**LCGW 
S  =  CM»LX (0*  »W) 

0=S*(0.2*S+1. Ci*(0.003333*S+1.0) 

N=2400* (C. l*S+1.0i 

G=N/D 

PEaREA  HZ) 

IMsAIMAG(G) 

FHI=57.2*ATAN2(IM»REJ 
IF ( PH  I • GT • C • )  PHI  =  PHI-360 
M  flG=C  A B  S  ( G  ) 

M0B=20.*ALCG1C(MAG) 

GREAL=L  !MlT(-£.t  5.,  RE] 

G I MAG=L IMIT(-5.»  5.,  IM] 

PH  1 1 30=  PH  I  + 18  C 

SAMPLE 

CALL  DRWG(l,l,LOGWtMCB) 

CALL  CP.WG(2»l.LOaWtPHJ 
CALL  DRWG< 1  ,2  ♦LOGW.O.O) 

CALL  0RWG< 2»2.L06h»-l80.) 

CALL  CRMG(3,ltGREAL1GIMAG) 

CALL  CRWG(4»1 »PHI » MDBi 
CALL  CP ’aG  (  2  » 2  tGPEALtC.O) 

CALL  CRtoG(4»2»L0GlA»Q«Q) 

TERMINAL 

CALL 

END 
STOP 


ENCPW(NPlCT) 
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TITLE  F'  =  3’J“‘::v  RESC1NSE 
TITLE  90  CE 

D  COMPLEX  S*NvCtGtCliC2 

0  CCMM0N/CAREA/S,N,LfG,ul,C2 

INTGER  N PLOT 
CONST  NPLnT=I 

CONTRL  FINTIM  =  1.C ,0  6LT  =  C.0005  ,DELS=0.0005 
PRINT  4.0E-3,W,PH,RE,  1 Y »  MOB  t  MAG 
D YNAM I C 

*”THE”SYSTEM  V ITH  THE  RESCNANcl  PEAKS. 

* - - - — - - - - -  —  - - - — - - - - 

PW=RAMP (C.C) 

LCGW  =  3  •  C  +  RIa 
10.**LCGVi 
S=CM°LX(0. »W> 

C1*S*(  C  .2*§4l  .CI*(C.G03333*S+1.0I 

C2  =  ( S*(0.2  5E-£*S  +  C.5E-5)  ♦i.0)*(S*<0.4E-7*S+0.632E  » 

-a  )*i  .o » 

N=24CC* ( C. 1*S*1.C ) 

r*  *  r  1  «.  j* 

g=n/6' ‘ 

PE  =RE  A  L ( G ) 

IM=AI MAG ( G  ) 

FHI=57.3*ATAN2(IM. RE  ) 

Ic ( OH  I .G7.C.  )  PH  I =  PH  I -3  6  C 
YAG=CABS(G) 

PCe=20.*AlCGlC(MAG> 

*  GREAL=L I  MI  T (- f • t  5.t  RE) 

*  GIMAG=LIMIT<-f .»  5.»  IM> 

FHIiaOPHI  +  l8C 

SAMPLE 

CALL  CRhG<lfI*LOGW,MCB) 

CALL  GRwG(2fltLGGw»PHI) 

CALL  CRV»G<  l,2iL0GW,0.C) 

CALL  CRWG(2,2tL0GW,-180.) 

CALL  CR!»G(3,lfP6.n'l 
CALL  CRfcG(4,l,PHl,MDB) 

CALL  0RUG(3»2»RE«Q.C) 

CALL  DRVG(4,2 tLOGWtO.O) 

TERMINAL 

CALL  ENCFMNPLCT) 

ENO 

STO° 


««* 


Till.'  FR  EQ!JCN CY  RESPONSE 
TITLE  BODE 

D  COMPLEX  S»N,0fG.Dlf02tC3 

D  COMMON/ CAREA/S,N,  C.G.D1.C2.03 

I NTGER  N PL.0 T 
CONST  N  PLOT  =  1 

CCNTRL  F INT I M  =  1  •  C  ♦  CELT  =  1 .OE-2 ,DELS=1 . CE-2 
PRINT  4.0E-3tW,PHl  ,RE » IM , MD8  ,MAG 
DYNAMIC 

PW=RAMP (C.C) 

LCGW=  3 •  C+RW 
W=10.**LCGW 
S=CMPLX (C. .VO 

C1  =  S* ( C .2*  S  +  l .0)*(0.003333*S«-1.0) 

C2S( S*<  C  .4E-7*S+0 •6324E-8J  +  1.0 )*( S*( 0 •25E-6*S  + 
C.5E-5I+1.C) 


*  - - - - - - - - - - - - apt 

*  COMPLEX  POLES  AT  FRECENCY  2800  HZ  InITH  D. RATIO  0.1  * 

*— — - - - - - - * 

C3  =  (S*(  127.55  1E-9*S+71.428E-6)*1.0) 

C=C1*02*C3 


IMAGINARY  ZERCS  AT  FRECUENCY  1900  HZ 


N=2300.0*<Q.1*S  +  1.0  )*(S*(27  7.008E-9*S+O.0K1.3) 

G  =  N/0 

Rc=P=AL(G) 

IM  =  AIMAG(G  ) 

PHI=57.3*ATAN2(IM,RE) 

IF(PHI.GT.C.I  FHI=PHI-36C 
MAG=CABS(G) 

M09=2C.*ALCG1C(MAG) 

GR EAL  =L  IMI  T< -  2  .  ♦  2.»  RE) 

GIMAG«LIVIT(-2.,  2  • .  IM) 

°HI180=cHl+18C 

SAMPLE 

CALL  OR V*G  ( 1 . 1 .LQGW.MCB  ) 

CALL  CRWG(2,1  .LOGW.PHI  ) 

CALL  CRWG< 1.2.L0GW.0.0) 

CALL  CRWG(2,2.LCGW,-18C.) 

TERMINAL 

CALL  ENCFW(NPLCT) 

END 


♦  TIME  RESPONSE  FOR  A  SYSTEM  CCHPtNSATEC  USING 
♦IMAGINARY  ZEROS  a  CCMPLEX  PCLES 
E°RCR=I  N-CUT 
IN  =STEP  ( C  •  C ) 

X8*ERRCP>»23C0  .C 

X7=REALPL<  C.0tC.3  333E-2»X8» 

X6=LE0LAG<0.1.C.2»X7) 

X5*CMPXPL(C .0 ,0.C» 5.GE-3»2000.0»  X6) 
X4=4.0E+6*X5 

X3=CMPXPL(0  .0,0.0 »  15.8  IE-6 » 5C00 . C» X4 > 

X2  *25 . C  E+6*X3 
STORAGE  CEN ( 3  )  • NUM ( 3 ) 

TAELE  CEN ( 1-3  J  =  8  3.3  33E-6,173 .611E-9, 1.0, .  .. 
NUM{ 1-3) «0.0* 277.CC8E-9tl.0 
X1=TRANSF(  2iDEN»2»MJMtX2J 
CUT=INTGPL(C.C»X1) 

TIMER  FINTIM=C.5tCLTCEL=C.005 
frtplt  cut 

LABEL  COMPENSATION  WITH  I.  ZEROS  C.  POLES 
LABEL  PCLES  AT  ?400  ZEROS  AT  1900  O.R.=C.l 
OUTPUT  TIME.OLT 
LABEL  S.  CONST  AIM  DCU  LAK  IS 

LABEL  POLES  AT  24G0  ZEROS  AT  1900  D.R.=0.1 
PAGE  X Y PLOT 
ENC 
STOP 
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TITLE  FREQUENCY  RESPONSE 
TITLE  COMPENSATION  WITH  COMPLEX  POLES 
0  COMPLEX  SfN,D,G,01,D2f02 

0  COMMON/ CARE A/ S,N,C,G,D1,C2, 03 

INTG5R  NPLOT 
CONST  NPLCT® 1 

CCNTRL  F INTIM® 1«  2  »OE LT =0.0006 1  ,D£LS=0 • C006 1 

PRINT  4.0E-3,W,PHI  ,RE,  IM.MDB.MAG 

DYNAMIC 

Rw=RAMP(C.0> 

LCGW=2 • 3+RW 
V=13.**LCGW 
S*CM°L  X  ( C. ,W) 

Cl=S*(0.2*S+l.C)*(0.003323*S+l.C> 

D2  =  ( S*<  C.4E-7*S+0.6  324E-a)  +  l.G)*(S*<  0.25E-6*S+, ..  . 
C.55-5J+1.C 


*  COMPLEX  POLES  AT  FRECENCY  19C0  HZ  WITH  0. RATIO  0.1 


* 

* 

SAMPLE 


C3*( S*<  277.008E-9*S+10.526E-5  )+l  .0) 
C=ri^P2*P3 

N=231C.Ci( C.1*S«-1.0  ) 

G  =  N/9 

RE  =  RE  AL  ( G ) 

I M=A I  MAG ( G ) 

PHI  =  57.3*ATAN2(IM,RE  ) 

IP(?HI.GT .  C  .  )  PHI =PHI -36C 
MAG=C AB  S ( G I 
MD8=2C.*ALCG1C<MAG) 

GREAL=LIMIT{-2.,  2.,  RE) 

GIMAG  =  LlVIT(-2«t  2 . »  IM) 
phub:=fhi  +  i8C 


CALL  CRWG( 1 ,1 •LQGto.MCB) 
CALL  CRWG(2,1 ,LCGW,PHI ) 
CALL  DRWG(1 ,2 ,LOGW,Q.O) 
CALL  CRWG ( 2 ,2  » LOG a  »-16C  •  ) 
CALL  OPWGt 2,1, RE,  IM) 

CALL  CRWC-(2,2»RE»C.O) 


TERMINAL 

CALL  ENCPW(NPLCT) 


END 


STOP 


TITLE  STEP  RESPONSE  CF  ThE  CCMPENS AT cD  SYSTEM  *ITH  u.PL'LcS 
INTGER  N  PLOT 
CONST  N  PLOT=  1 

CCNST  K 1=25. CE+6,P 1=31. CE-6 , P2=500G.C 
CONST  K 2=4. CE+6, P3  = 5. CE-3»P4= 2000.0 

CONST  K2=4.2C25E+6,P5=C.3333E-2,P6=0.  1 , P7  =  C. 2 , P8  =  0 . C5 
CCNST  P9=23CC.C»K4=2300  .0 
I NT EG  R  KSC  X 
DERIVATIVE 

E=R-Y 

R  =S TE P ( C  •  C  ) 

Q=CMPXFL(C  .0  »0.0.P1.P2»K1*EJ 
P=CMPXPL<0.C,C.0,P3,P4,K2*GJ 
M=REALPL(C.0,P5tP) 

L=LECLAG(C.0,P6,P7tM) 

YCOT=CVPXPL( C.C,C.C»P8,F9,L*K3J 
Y=INTGRL(C.0,YDCT*K4) 

CCNTPL  P INTIM=  0.7,CELT=2.CE-4,DELS=4.CE-4 
PRIN~  C.Cl»E,C.YDCT,L,Y 

SAMPLE 

CALL  DRWGI1»1*TIPE*YI 
TERMINAL 

CALL  ENCPW(NPLCT) 

END 

STD0 


TITLE  FO'Q'JENCV  RESPONSE 
TITLE  8C  CE 

D  COMPLEX  S,N»C,G»C1»D2»C3 

0  CCMMON/CAREA/S  ,N,C,G, 01,02 ,D3 

INTGER  CUR » N PICT 
CONST  CUR=l ,NPL0T=1 
PARAM  Kl=?06.6llE-S,K2=S0.9CSE-6 
CCNTRL  FINTIM*l.C,CELS*6.25E-4 
PRINT  4.0E-3  ,W , PHI  ,FE » iM *  MOB ,MAG 
DYNAMIC 

PW=RAMP (C.O) 

L0GW=3. C+RW 
W=10.**LCGW 
S=CMPLX (C. ,W) 

C1=S*(C.2*$+1 .C)*<C.003323*S+1.C) 
C2=(S*(  C.4E-7*$4-C.6324E-8)+1.0)*(S*<  C. 


25E-6*S+0.5£-5)+ 


COMPLEX  PCLES  AT  FRECENCY  2200, 220C.25CC  HZ  WITH  0. RATIO 

C3=(S*<K1*$+K2)+1.C) 

C=C1*C2*C3 

W  ^  (  w  •  X  f  O  ^  X  •  1*  ) 

G=N/D 

RE=REAL (G) 

IY=AIMAG<G) 

PHI=57.2*ATAN2(IM,RE) 

IF(BHI.GT.C.)  FH-PHI-26C 
MAG*CA8S (G ) 

M0B=2C.*ALCG1C<MAG) 

C-REAL  =  LlMIT(-2.,  2.,  RE) 

G I  MAG  =L I WI T  < -2  • »  2.,  IM) 

FHliaO=PHl+18C 

SAMPLE 

CALL  DRWG(l,CLRf LCGW,MC8I 
CALL  CRWG(2tCLP,LCGW,PHI  ) 

CALL  DRWG(1,4,LOGW,0.0) 

CALL  DR*»G(2,4,L0GW,-18C.) 

TERMINAL 

IF (CUR. EC. 3)  CALL  ENCRW (NPLOT ) 

CUR=C  UP ♦ 1 

ENO 

PARAM  Kl*ia<9.C35E-9,«2  =  8fc.S56E-6 
END 

PARAM  Kl=16C.CCCE-9,K2*8C.000E-,6 

END 

STOP 


79 


LIST  OF  REFERENCES 


John  J.  D  1 ftzzo/Const  ar.tine  a.,  Houpis  Linear  Control 
System  ft r. aly sis  ani  Design  £114  Modern 

McGraw-Hill,  1982. 


Thaler,  G.  J.,  Das ian  of  ? aaib ack  Systems,  Dovden, 
Hutchinson  5  Ross, la z  ,  1  97 3 . 


Ogata  K.,  Modern  Control  Engineering,  Prentice  Hall, 
Ir.c,  1970. 


Syn,  W.  M. ,  Turner,  N.  N.  aal  flyman,  D.  G. ,  Digital 
Simulation  Language  J  ser  Manual,  I3M  System/36 o7~19687 


Di  Stefano.,  Feadoack  and  Control  Systems,  Schaum 
McGraw-Hill,  1  9807 


80 


INITIAL  DISTRIBUTION  LIST 


No.  Copi 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  California  93  940 

Department  Chairman,  Code  5  2  1 

Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93  940 

Professor  3.  J.  Thaler,  Code  62Ir  5 

Department  of  Electrical  Enginaering 
Naval  Postgraduate  School 
Monterey,  California  33940 

Professor  A.  Gerba,  Jr.  ,  Code  S23z  1 

Department  of  Electrical  Engineering 
Naval  Postgraduate  School 
Monterey,  California  93  940 

K.  Constando ulakis  1 

Stratigcu  Rogakou  49 
Polidrossor.  Amarousioi 
Athens,  Greece 

TNFG-UN  Raul  E.  Samaniego  1 

P.  0.  Box  95  3 
Guayaquil  -  Ecuador 
South  America. 

Sozon  A.  Constandoulaki s  6 

Andoniou  Paraskeva  5 

Amaroussion 

Athens-Greece 

Emmanuel  Hor iano poulos  3 

Naval  Postgraduate  School, 5 .H. C. 1 363 
Monrerey  CA- 93940 

Helenic  General  Naval  Staff  2 

Education  Department 
Stratopedon  Papagou,Holargos 
Athens  Greece 


81 


